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Abstract 

The thermal propagator of a heavy muon propagating in a hot QED plasma 
is examined within the Bloch-Nordsieck approximation, which is valid in the 
infrared region. It is shown that the muon damping rate is finite, in contrast 
to the lower-order calculation with hard thermal loop resummations taken 
into account. 
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I. INTRODUCTION 



Because of much interest in the early Universe as well as in the quark-gluon plasma 
to be produced in heavy-ion collisions, many authors have studied gauge theories at high 
temperature. Thermal propagators of bosons in real-time thermal field theory carry the 
Boltzmann factor l/(e^/"^ ~ !)• For a massless boson like gauge boson, this factor behaves 
as ~ T/k at the infrared region, k <^ T, which leads to stronger infrared divergences 
in thermal amplitudes than in vacuum theory. In vacuum theory, it is well known 
that all the physical quantities are free from infrared singularities, provided one sums over 
all degenerate final states (Bloch Nordsieck mechanism). An important exception to this 
general statement is the so-called Coulomb singularity. For thermal reaction rates, due to the 
stronger infrared singularities mentioned above, it is not known if they do not possess infrared 
singularities, except again for Coulomb singularities. Up until now, a proof is available of 
absence of leading infrared singularities in a generic thermal reaction rate |^ and of all 
infrared singularities in a bremsstrahlung process within some specific approximation [Q. 

It has been established [|,^|14| that, to obtain a consistent perturbation expansion, 
resummations of hard-thermal-loop (HTL) contributions are necessary. According to this 
HTL resummation scheme, infrared behaviors of thermal amplitudes are softened or screened 
compared with those in naive perturbative calculation, this softening renders some otherwise 
divergent physical quantities finite [0)11 • There are, however, some other physical quantities 
which are still infrared divergent to leading order in HTL resummation scheme. A well- 
known example is the quark (muon) damping rate in a hot QCD (QED) plasma 0. As 
a matter of fact, the quark (muon) damping rate diverges due to the Coulomb singularity. 
In naive perturbation calculation, the divergence is of power type. In HTL resummation 
scheme, as mentioned above, the infrared behaviour is softened and the infrared divergence 
becomes to be of logarithmic type. It is expected, at least for quark case, that the infrared 
singularity is screened when still higher order resummation is performed. As far as the 
infrared behaviour of a spinor field in QED is concerned, Bloch-Nordsieck approximation is 
a good approximation [0,11|. In this paper, the thermal two-point function of a heavy muon 



in a QED plasma is computed within this approximation. Various properties of the result 
are discussed, among those is an especially important conclusion: The decaying behavior of 
a heavy muon does not suffer from divergence. 



II. THERMAL TWO POINT FUNCTION OF A HEAVY MUON 

By the term "heavy muon" we mean so heavy a muon that it is not thermalized, i.e., 
g-s/T ^ rjj^g j^^^ QED plasma means that the temperature T of the plasma is much 
higher than the electron mass m, T ^ m. As discussed in Sec. I, the infrared behavior of the 
theory is responsible for the divergent damping rate. Thus, we employ the Bloch-Nordsieck 
approximation to the muon sector of the QED Lagrangian. The Lagrangian density is 
obtained [10|JTT|] by substituting a constant vector (u^ = 1) for 7^ in the muon sector of 
the standard QED Lagrangian density: 

£total _ £^ _|_ £^ _|_ £^ ^ (2.1) 
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= ^{iu -d + eu- A- M)^ , 
= ■ d + e'-y ■ A — m)ip , 

F^^ = d^A^ - dy 



A 



(2.2) 
(2.3) 

(2.4) 
(2.5) 



where M,m, and A stand, respectively, for heavy muon field, electron field , muon 
mass, electron mass and photon field. It is known that, in the infrared region, the spinor 
structure of the heavy muon does not play any significant roles |1T0|JTT|1 . Then the above 



Lagrangian density well approximates the exact QED Lagrangian in the infrared region. 
To obtain the thermal two-point function of the heavy muon, we start from the Euclidean 



or imaginary-time formalism of thermal field theory, initiated by Matsubara [T^, in which a 
time arrow fiows from zero to — i/T in a complex time plane. In this formalism, the energy is 
discrete; = 27mT {n = 0, ±1, ±2, ■ ■ ■) for a boson field and = 7r(2n + 1)T for a fermion 
field. At the final stage, we analytically continue the two-point function, thus obtained, 
to real time (or real energy). A four vector a'^ = (0,04) in the imaginary-time formalism 
is obtained from the counter part (oq, a) in Minkowski space by an analytic continuation, 
ao — > —ia^. Similar continuation is made for tensors. 
The generating functional is introduced as usual. 



V^V^VipVipVA exp S . 

^^xfc'"'''^ - \{d- Af 
V 2A 



Z 

s 



where the covariant gauge is employed. In Eq. (p.7|) , rj, ■ ■ ■ are sources being 
to the fields with which they couple. After integrating over \E', \E', ip, and ip in Eq. 
two point function of the heavy muon is obtained as. 



(2.6) 



(2.7) 

conjugate 
(pj), the 
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ri,-=0 



Z 6ri{x) 6r]{y) 

^ J VAdet{iuE ■ d - euE ■ A - M) 

X det(z7 ■ 9 — 67 ■ A — m)G{x — y; eA, M) 



X exp 



2A 



id ■ AY 



(2.8) 



where 



{iuE ■ d — euE ■ A — M)G{x — y; eA, M) = 6Eix — y) 



(2.9) 



5E{x-y)^T Y: I 



ki=odd 

h = vr(2n + l)r 



3 exp[ik ■ {x - y)] , 
(n = 0,±l,±2,---) 



(2.10) 
(2.11) 
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In Eq. ( |2.10| ) , "fc4 = odd" indicates that the summation runs over n as defined in Eq. ( p. Ill) . 

Since the heavy muon is not thermahzed, the muon determinant has been reduced to 
di.ei{iuE-d — M) as in vacuum theory |Ty]. Using the formula det a = exp(Tr In a) and Furry 's 
theorem [Ol, we have 



det(i7 ■ d — e'j ■ A — m) = det(i7 ■ d — m) exp 

L Z J 

x(7-A)Go(z2-^i;m) + 0(e^) 

f d^k exp[ik ■ {x - y)] 

Go{x-y;m)- ^ 



6 f 

—Tt J d'^zid^Z2{'j ■ A)Go{zi - Z2; m) 

(2.12) 

(2.13) 



, ^ J (27r)3 + m 

Here Go{x — y]m) denotes a bare electron propagator. As seen from Eq. ( |2.8|) with Eq. 
( p.l2| ), in the heavy-muon two-point function, the contributions of electrons enter through 
thermal electron loops. Here we employ the HTL resummation scheme for the electron 
sector: To leading order, among such electron loops, only HTLs that yield corrections to 
soft-photon external lines are important. In hot QED considered here, the O(e^) term in 
Eq. ( |2.12| ) vanishes [1^1 to leading order in the HTL resummation scheme. We then ignore 
it in the following. 



Following we introduce the integral representation, 



^00 

G(x — y; eA, M) = —i drexp 

Jo 



iriuE • d 



-CUE ■ A- M + iec) 



Se{x - y) 



where ec is a "convergence" factor. Defining the function U{t) as. 



U (r) = exp 



iriiuE ■ d — euE ■ A — M + itc 



Ssix-y) 



G(x-y]eA,M) = -% / drU{^ 

Jo 



we have 



dU 

— i—- = Uue ■ d — euE ■ A — M + iec)U(T) , 
or 

U{0)=6E{x-y) . 

Now we assume the following form for the solution to Eq. (|2 . 1 7| ) with Eq. (|2.1S 



(2.14) 



(2.15) 



(2.16) 



(2.17) 
(2.18) 



Uir) =T j 



d?q 

q4=odd ' 

, (2vr)3 



exp 



iK{x — y\A, r) 



+iq ■ {x — y) — ir^M + ue ■ q — i^c 
Eq. ( pT7| ) with Eq. ( l2J9| ) leads to 



(2.19) 
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dK dK , , , 

- — = UE--^ + euE-A, (2.20) 

K{t = 0) = 0, (2.21) 
where ue = {u^u^). It is straightforward to solve this equation in the form: 

—- exp[ik ■ (x - y)]A{k) ■ ue rfr'expHw^; ■ kr'] , (2.22) 



k4,=even ' 



^4 = 2TTnT (n = 0, ±1, ±2, ■ ■ ■) 



(2.23) 



where "^4 = even" indicates that the summation runs over n as defined in Eq. ( p.23|) . 

Using all the formulae displayed above, and carrying out the integration by part, we 
obtain for in Eq. dJ): 



iq ■ {x — y) — iT{M + ue ■ q — i^c) 



qi=odd ' 



6 f 

— — j d'^zid'^Z2TT (7 ■ A)G{zi — Z2;e = 0, m)(7 ■ A)G{z2 — Z\\e = 0, m) 

—ieT / A{k) ■ Ue expHk ■ {x — y)] [ dr' exp[—iuE ■ kr'] 

1 _ J [2t[Y Jo 



Going to the momentum space, 

Qe = j d^{x - y) exp[-zp • (x - y)]GE 



(2.24) 



(2.25) 



we have. 



Qe = i I T^A I drexp 

dn 



-ir^M — Ue ■ p — i€c) 



E 

l4=even 



/d /"^ 
— exp[iZ ■ uet\A{1) ■ Ue dr' exp[-z/ • uet'] 



(27r)3 

dH 
(2vr) 



where 11^'^ is the thermal vacuum polarization tensor of the photon: 

d^k „ „ 1 .,1 



(2.26) 



7I 



k4=odd 



(27r)3 ' ^.{k-l) + M' -i-k + M 



p 



(2.27) 
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where 5^'^ = diag (1, 1, 1, 0) and = (1,0). In Eq. Q^ ), Ut (Ul) is the transverse 
(longitudinal) component of 11^'^ (See Appendix A). 

The exponent in Eq. ( |2.26| ) is quadratic with respect to A, and we can perform the 
functional integration over A: 



dr exp 



a 



-iT{M + UE-p- iec) + —^D{t, ue) 

ZTT^ 



(2.28) 



D(r,Mi5) = -(27r^)3-^T ^ f d^'l —^-—{l - cos{l ■ uet)] 



X 



+A 



I ■Ue\ 



+ 



2 (I-Ue)^ ^2 , (^■^) 
Ue 77^ U + 



P 



(2.29) 



where a is the fine structure constant. Here, in order to regulate the ultraviolet divergence, 
we have employed the dimensional regularization, D = 3 — e with e > 0. (A dimensionful 
parameter is introduced, as usual, so that e remains dimensionless.) In the HTL approxi- 
mation, which we employ for the electron sector, the explicit form of 11^ and Ul are known 

mm-- 



1 2^4 . , 

--nisc — \ 1 

2 P 



il 



l + -^\-LL{th,\l\) 



21 



Ua 



1--J|^(^^4,|/|)| , 



(2.30) 
(2.31) 



where 



L{ih, \l\) = In 



eT 

71' 



ih - \l\ 



(2.32) 
(2.33) 



nisc in Eq. (|2.32|) is thermal or Debye mass of the thermal photon. In obtaining Eqs. (p.29| ) 
- (|2.33|) , the electron mass m has been ignored (cf. the first paragraph). 

Let us turn back to D{t,ue) in Eq. (|2.29| ). After a straightforward but tedious calcula- 
tion (cf. Appendix A), we have 



U4 



D{t,ue) = -2tt'Tt^ - + ln 



\u\ 



27Tfi 

VT^rris 



1 



g-msc|w|T 



\U\Tmsc 



+ Ei{-\u\Tmsc) 



-2'n^TT\u\ 



1 7 , 2Tin,\u\T 



(A - 1)-k'^Tt\u\ 



+2n'TT 



2^ 



\u\ 



7 



27r 



— — In V TT H h 7r7 — 27r In 



Ati^PT 27rM4 
2nfj, \u\ 



tan 

\u\ 
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In h In 



TTUe 



\U\[\U\ — 2M4j 



27111 

ln(l - e 



+ 47r^Tr- 



ax— 



1 



^2-kT\u\tx _ ]_ 



dx- 



-2TTTT{\u\-iUi) 



1 



) + C.C. 



1 







+ C.C. 



+vr(A-l) 
— 2|-u|r7rT 



i _ In ^ + ^ + (-^ + ln47rT) - 2 ln(l - e'^^^l^^) 



+ C.C. 



g27rT|«|r _ 1 ' 2 

7rTr(|n| — 



1 y. iln(l - g-27rTr(|«hm4)j 



g27rTr(|u|-m4)-l 
1 



+ C.C. 



g27rTr(|«|-m4) ]^ 



+ C.C. 



+3m47rTr / dx 
Jo 



27rTr(|ii|— m4a::) 



(2.34) 



i\—iu4x) ]^^2 

where 7 ~ 0.57- • ■ and Ei{—\u\mscT) are, respectively, the Euler constant and the exponen- 



tial integral function and ue = yU4^ + 1^1 • The details are given in Appendix A. 

An analytic continuation of in Eq. ( |2.28|) with Eq. (|2.34|) back to Minkowski space, 
which is obtained through ^4 — ^ i{pQ + ze), — > i{uQ + ie), yields the retarded function. 
We normalize it through MS scheme ||^ , which amounts to choose the renormalization factor 
Z2 as 



„ 3-A 1 

Z2 = e"—- 



Then we obtain, 

QR = iJ^ c/rexp |-ir(M - M -p - kc) + ^r"/?! 
where u ■ p = u^p^ is the scalar product in Minkowski space. 



(2.35) 



(2.36) 



III. PHYSICAL CONTENTS OF THE TWO POINT FUNCTION 

We study the behaviors of in Eq. (p.36|) in various limits of parameters involved. 
We divide the region < r < 00 into several regions, which are discriminated by the two 
parameters ta = tb = (see Appendix A). 

a) r ^ 

b) < r < tb 
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+27r- 



2 I2 II , 



vr— In 

\u\ 



Uo - 


\u\ 


Uo + 


\u\ 



Uo 



\u\ 



(3.2) 



c) < r 



Vr ^ -271^Tt\u\ J In tT- — + 7-1-^ 

2 MM 



In 



T 



-2vr^T^ 



+ 27r- 



i[2 



lulm^c 2 



A-1 



In 

2 \2 /i / 



TT— In 

m 



Uo - 




Mo + 


\u\ 



2 Mo 



(3.3) 



In the above equations, use has been made of = 1. 

When the muon is far from the mass shell \M — p ■ u\ ^ M, in Eq. (|2.36|) , the region 
T ^ gives the dominant contribution: 



Gr 



^r(i + o 

M — u ■ p 
3- A 

5 



i{M — u ■ p) 



a 



27r 



a 



A = exp — < 1 

vr 



3- A 
4 



(7 — Invr — ivr) 



(3.4) 
(3.5) 
(3.6) 



where F is the gamma function. When A > 3, < in Eq. ( p.4| ). This is just as in vacuum 
theory |jlO|, which we do not discuss any further. When A = 3, in Eq. (|3.4| ) behaves as 
the bare propagator. 

We are now in a position to discuss the decaying behavior of the heavy muon. For 
this purpose, we should study Qr near on the muon mass shell \M — p ■ u\ ~ 0. The 
time dependence of the heavy-muon two-point function is obtained through the Fourier 
transformation; 



Ci? = y ^ exp[-ipot] 



— exp[—iEt\ exp 



a 

2^ 



■Vnir = Mt/E, u) 



(3.7) 



In obtaining Eq. p.7| ), use has been made of the on-shell approximation {M + u- p)/uo — E. 
From Eq. ( p.7|) and Eqs. ( |3.1| )- (|37^ ), it is straightforward to obtain the time dependence of 
Gr- 

a) T <tiTA 



e' [v^e) 



A 



(7 — in) 



(3.8) 



b) TA < r < TB 



8 



E 



E 



E 



\u\ 




Uq 



/i / 2\u\ 



In 



Uq - 


\u\ 


Uo + 


\u\ 



— in- 



2\u\ 



(3.9) 



c) < r 

m 



^ iEt+F{t) 

E 

MTt,^, f, MTt in 

— a \U\ < In + 7 — 1 



(3.10) 



Mo 



T 



Mm, 




Uo - 


\u\ 


Uo + 


\u\ 




We discuss two cases separately. It is to be noted that |-z2| = = ^^^J^i— 
the external muon velocity: 

• <v <1 

In the region a) above, the behavior of Gr is gauge-parameter dependent, and we can- 

not draw any physically sensible conclusion. In the region b), the factor ( ^ j 

in Gr, Eq. (|3^) , dominates over others. Then, the leading part of Gr is gauge- 
parameter independent and increases as t increases. In the region c) the dominant 
part of Gr is gauge-parameter independent, which first increases and then decreases. 
The turning point in t is approximately given by ■^{ ^'^\^\ Y^^ ■ 



In the region a) and b) the situation is the same as the first sentence in the above case. 
In the region c) Gr decreases monotonously. 

Let us discuss the large t region c) more elaborately. The leading contribution comes 
from logarithmic parts of F[t) in Eq. ( p.llD (cf. Eq. ( |3.3| )). It is to be noted that the 
gauge-parameter dependent part in Eq. ( |3.11D may be neglected and Eq. ( p.ll| ) is gauge 
independent: 



F[t) ~ —a — —\u\ 
E 



In 



MTt 
E 



Uo 



\u\ 



In- 



T 



m.rU 



(3.12) 



A striking feature is that Eq. (|3.12|) is divergence free or finite. The damping rate 7^ is 



usually defined through the behavior of Gr_ at large t; Gr 



~ e 



it 



00). Ware it 



not for the factor Int in Eq. ( |3.12| ), the damping rate would be identified from Eq. ( |3.10| ) 
with Eq. ( p.l2| ). Eq. ( |3.12|) indicates that, at large t, the heavy muon in a QED plasma 
decays much faster than the constant damping rate. In any case, the quantity, F{t), in Eq. 
( p.l2| ) governs the damping behavior of the heavy muon. The behavior of F{t) in Eq. (|3.12|) 
is as follows: 
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• < < 1 



F{t) decreases after t = td{= -^ (^^ ^Y^^"^)- This means that after t^, the heavy 
muon begins to decay. 



• t> ~ 1 

In this case, 



Fit) ~ -a 



MTt , Mm,rt 



(3.13) 



This is always negative and the heavy muon decays. 



IV. SUMMARY AND DISCUSSION 



In this paper we have derived and discussed the two point function of a heavy muon 
passing through a QED plasma within the Bloch Nordsieck approximation. The assumption 
of heavy muon means that the muon line is hard. The electron sector only comes in as 
thermal radiative corrections to photon lines. Among those, we only take into account the 
HTL corrections to soft photon lines. HTL corrections to vertices are not necessary. 

It is well known that, to leading order of HTL resummation scheme, the damping rate 
of a quark (muon) diverges due to the Coulomb-type infrared singularity even to the lowest 
non-trivial order of resummed perturbation theory. It is because only the electric sector of 
soft gluon (photon) line is screened by the thermal mass. As far as QCD is concerned, the 
infrared singularity is screened by the magnetic mass which is expected to arise through still 
higher order resummations. But in QED, no magnetic mass is induced. 

In this paper, within the Bloch Nordsieck approximation, it is shown that the "damping 
rate" of a heavy muon in a hot QED plasma is finite. More precisely, the heavy muon decays 
as at large t. 

Finally it is worth mentioning the case where the HTL corrections to the soft photon 
line are ignored. (Explicit computation is given in Appendix B). For \M — p ■ u\ ^ M, the 
gauge parameter dependent parts of Qr (cf. Eq. ( p.36| )) is the same as Eq. (|3.5|) . Also for 
A = 3, ^i? reduces essentially to the bare propagator. In the nearly on the mass shell case 
the time dependence of Gr is quite different from that of Eq. (|3.10|) with Eq. (|3.12|) : 



the anti-damping behavior! From this fact, we learn that the HTL resummation for soft 
photon lines plays an crucial role for obtaining a physically sensible damping behavior. 
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APPENDIX A: CALCULATION OF D{t, Ue) 



According to the HTL resummation scheme, when the external photon hnes are soft 
\l\ ~ O(eT), the HTL corrections are necessary. Since 1^ is 27inT {n = 0, ±1, ±2, ■ • •), 1^ 

with n 7^ is no longer in the soft region. Then we should take care of the n = mode only. 
The explicit expression for IIt,l in Eq- (|2.27|) , within the HTL approximation, are given 

in [^,15|. From these expression, we have 

Hr(n = 0)=0, ULin = 0) = mj = ^ . (Al) 
We now decompose V{t,ue) in Eq. ( p. 291) into two parts: 

V{t, Ue) = -T{27Ti2y ^ J d^'l V{t, Ue; h, 

li=even 

= -T{2'K^iy J d^l P(r, Ue; h, 

-T{2nfiY Y.' Jd''lV{T,UE;U,l) , (A2) 



li=even 



where the prime on the summation symbol in Eq. (^^) means to take summation over /4 
with non-zero modes. We further decompose Eq. (|A2|) as 



- T{27rfiy j d^'l V{t, ue; k, I) = '^rn , (A3) 

•' m=l 

-T(27r^)^ Y' J d^^l V{T,UE;k,l) =j2l^'m, (A4) 



where 



l4=even m=l 



{l-uf P + mj 



^ 1 — cos(/ ■ ut) u^ 



V, ^ -T(27r/i)^ / d^l \ ' ^ , (A6) 

(/ ■ uy 



Y). 1 — COs(/ ■ Ut) I I ■ U 



v,^-nx-m.^Yjd"i . (A7) 



l4=even 



{I ■ Ue) 



l4=even 

2 



(/ ■ Ue) 
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It is to be noted that the integrand of Vm, {jn = 1, 2, 3) and V'm, ijn = 1, 2) are regular at 
r= 0. 

We first evaluate V^, [tti = 1,2,3). Carried out the integration over the polar angle 
cos~^(||i||), Eq. (^) becomes 

V, = -T(27r/i)^j|^27ri-^/2r(e/2) dl {f + ml)-^ I + c.c.j . (AlO) 

In order to perform the integration over /, it is convenient to analytically continue the 
integrand into a complex I plane. By deforming the integration contour (as shown in Fig. 
IID, we obtain. 



Pi = -T{27rfiY^27r^-'/^T{e/2 



2 



\U 



71t\u\ _^ i f°° VP - nisc^) 



SC j J 



dlS^ e-^f ^(1 - e-l">') + c.c. 



_^ Urn / dl ^^' ~ "^'/^'^" {l - e-l">') + C.C. 



rmsc 


dl 


Jp 






(I 


\u\ 













27ru 7 1 - e-l"l™- ^ , , ^, 

+ Ei(-\u\Tms 



+ In 



^/7^msc 2 \u\rmsc 
2-71 fi 7\ 



^fnvfisc 2 ) 

{|M|r(l — 7 — In \u\Tmsc) + Od-upr^m^c) (l-ulrmsc <^ 1) 
T . (All) 

— + 0(^^) (l^lrm., » 1) 

Through similar procedures, T>2 and T>^ may be evaluated: 



V, = -2n'T\u\r [- + \n + ^ - 1 (A12) 

V e ^/v^ 2 / 



Ir I - + ln ^L' + - - 1 

Vs = -(A - l)7r^T\u\T . (A13) 

The terms proportional to 1/e , which correspond to UV divergent piece at e = or D = 3, 
is to be eliminated later by renormalization. 

Let us evaluate V'm (m = 1,2) in Eq. (|A4|) . To evaluate both Eq. (|A8|) and (|A9|) , it is 
convenient to introduce 



The integral in Eq. ([A8|) is obtained from F(l4) with a = 0, while the integral in Eq. (|A9|) 
is obtained from dF{l^/da with with a = 0. Carrying out the integrations over the polar 
angle and using the technique of deforming the integration contour in the /-plane, we obtain 



F{U) = -T{2n^yn^-'/^T{e/2) 
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+ sin^ r dl {I'-ill + a)}- 



e/2 



1 



{u4i + i\u\iy 



^iU4l4T—\u\lT ") 

+ C.C - + C.C 

[U4I4 + t\u\l)'^ J 



(A15) 



As to the integration of the third term in last brackets in Eq. (|A15D , we first carry out the 
integration by parts and the integration by deforming the contour Ci to C2 C3 in Fig. 2. 
Then we shift the integration vahable to get 



F{k) = -T(2TiiiyTi^-'/^T{t/2) 



\u\ \ \ur 



-6/2 



a 



+ sm y(/4 + a) 



2 / dl 
h 



{f - 1) 



-e/2 



il\u\ + 



lU4 



2 + C.C 



J- roo g-l\u\T ^il^u^T-\u\T^l\-^u 
+ 7-^ / , dl ; V . + C.C. 



+7-^ / dl — , h C.C. 

1^1 -^0 -nl + a-il 



(A16) 



Using Eq. ( [A 161 ) with a = and carrying out the summation over we obtain 



l4=even 



(7=0 



2tx^Tt^ - ? - In ^/tx] + — + 7r7 - 27r In ' ^ - tan"^ 



u\ Ve 2 
2 / y 



.2vr2rr^ I In 



+ 7r7 — 27r In ■ 
e 27r/i 

, 1 



■U4 
\u\ 



TTUe 



+ In ^ + Atc'Tt 

2TTfl \u\ I \U\ Jl 



dx 



X 



^2ttT\u\tx _ I 



^, ^, , ln(l - e-2-^"(l^l-™4)) + C.C. 
u\{\u\ — W4) 

^ lo \u\ -iu4,x e2'^^^(l"l-»«43;0 _ 1 



\u\ 



+ C.C. 



(A17) 



l4=even 



V'2 in Eq. ( |A9D may be evaluated as: 

P'2 = -(A-1)_ E' [T(2vr/ir/rf^/^ 
and we obtain, 
V'2 = -{\-1)tx 



dadr"^ 



F{h) 



u=0 



(A18) 



1 



+ i + In ^ - 7 + (7 - ln4vrT) + 2 ln(l - e--^l->) + 2\u\rnT -^^^^^^^ _ ^ 
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-lln(l-e-2-^^(l"l-^"^))-c.c. 

1 



— c.c. 



/•I 1 
—?>iuATiTT I dx ^ ^ . — r h c.c. 

/•l g— 27rrT(|u|— m4x) 

+z2r27r2TV dx {\u\ - iu,x) ^^_^^Tr[\u\-iu,.) _ + C-^- 



(A19) 



Now let us study the behaviors of T)' ^ (m = 1, 2) in the following two limits (r^ is as in 
the text, ta = j^): 

a) r -C 

I,',^2.flH-I + l+l„5^U2.Vr^fi + 2-l + lnB^) , 

\e 2 2^/7!■ J \u\ \e 2 a/tt y 

(A20) 

^ -(A - l)vr + I + 1 + In + (A - l)n'T\u\T , (A21) 

and, 

b) r > 

P'. ^ 2. f 1 + 2 - m ii!^ - tan- + (i - ^ + In ^) , ( A22) 

\e 2 2TTjj, \u\ \u\ J \u\ \e 2 UE\jT^i J 



1 7 , 47r3/2T\ 
7+2 



P'2^-(A-l)7r(- + ^-ln^-— j . (A23) 



Substituting all the results obtained above into Eq. ( |A2|) , we deduce the following be- 
haviors of T){t^ Ue) in various limits (r^ is as in the text, Tb = 

a) r ^ Ta 

.(...).2.(i.I.l.l„!2^)-.(A-I,(i.2..^) , ,A24, 

b) TA < r < Ts 

P r, Mb ~ -2n^TT-^ In rTu^ + 7- l+27r - + --ln ^ tan"^ ^ 

|m| \e 2 27r/i |m| 

-(A - l)7r2T|^|r " " l)vr + | - In j , (A25) 

c) Tb < r 
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Ue 



V{t,ue) ^ -27r^Tr^ (In rTuij + 7 - 1) +2 



'1 7 , 47r3/2T M4 



\u\ 



TT I - H In ■ 

e 2 27r/i 



tan 

Ml \u\ 



-(A-1) 



TT^T M r - A - 1 TT - + - - In 

' ' ^ ^ Ve 2 27r/i y 



-27r'^T\u\T—^(—'j + 1 — In \u\Tmsc) + 27r^r 



(A26) 



The remaining tasks are to eliminate UV divergence through renormahzation with MS 
scheme and to perform the r integration (see text). 

APPENDIX B: NO HTL CORRECTION TO THE PHOTON PROPAGATOR 



We set Ut,l = in Eq. (12^2^) . In this case, Eq. ( ^^291) reduces to 



V{t, ue) 








mac=0 


msc=0 



=0 m=l 



m.c=0 



{I ■ uY 

^ + (A-l)(^"^"' 

P 1/2 



X < 



Here 



V 



1+2 



-T {27V i^y J d^l 



J) 1 — cos(/ ■ ut) Ue^ 



msc=0 " {I ■ U 

= -(A - l)T{27Tfiy I d^'l 



msc=0 

Similar procedures as in Appendix A leads to 

2 



e / jD, 1 "~ C0S(/ ■ ut) I ■ U 



^1+2 



msc=0 



msc=0 



2vrVT^ i + ^-l + ln 
u \e 2 



{I ■ m)2 

27r/i|n|r 



/2 



-(A - l)7r^T|M|r . 



(Bl) 



(B2) 



(B3) 
(B4) 



(B5) 
(B6) 



Substituting the above Eqs. (|B|) and (0) and Eqs. (|A20D -(|A23|) to Eq. (^), we have 
a) T <^ ta 



V{t,ue) 
b) r > ta 



msc=0 
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V{t,ue) 



_27r2Tr ^ (In T^et + 7 -l) + 27r(- . ___ 

\u\ \e 2 ZTT/i 



msc=0 



7 , 47r3/2T _iM4 

— Ill — -— tan -— 



\u\ 



\u\ 



(B8) 
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FIGURES 

FIG. 1. Integration contour in the complex / plane 
FIG. 2. The integration along the contour Ci may be evaluated along the contour 6*2 C3. 
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Fig. 1 
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